Mode locking of vortex matter driven through mesoscopic channels 
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We investigated the driven dynamics of vortices confined to mesoscopic flow channels by means 
of a dc-rf interference technique. The observed mode-locking steps in the IV-curves provide detailed 
information on how both the number of vortex rows and the lattice structure in each flow channel 
change with magnetic field. Minima in flow stress occur when an integer number of rows is moving 
coherently, while maxima appear when incoherent motion of mixed n and nil row configurations is 
predominant. Simulations show that the enhanced pinning at mismatch originates from quasi-static 
fault zones with misoriented edge dislocations induced by disorder in the channel edges. 
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During the last decades, vortex matter in superconduc- 
tors has been widely recognized as model system for in- 
vestigating driven, periodic media in various pinning po- 
tentials. One of the intriguing phenomena it can display 
is a non-monotonous behavior of the depinning current 
density J c with magnetic field. The best known exam- 
ple of such behavior is the so called peak effect occur- 
ring in weak pinning superconductors near 77 C 2 JjJ . Non- 
monotonous behavior of J c also occurs in superconduc- 
tors with (artificial) periodic pinning arrays due to the 
possibility of (in)commensurability between the pinning 
structures and the vortex arrays In both cases, the 
change in J c is accompanied by pronounced structural 
transitions in the static as well as the driven " lattice" . 
For periodic pinning, it was predicted theoretically 0] 
and confirmed by experiments Q] that defects due to in- 
commensurability can cause a drastic reduction of J c . 
For the case of the peak effect in presence of random 
pinning the consensus is that the opposite happens: here 
plastic deformations cause an increase in J c . However, 
the role of defects in the plastic flow regime above the 
pinning threshold remains largely unknown. 

A powerful tool to obtain time averaged, small scale 
information on the moving structure is to probe the dc 
current- voltage (7- V) characteristics in presence of super- 
imposed rf-currents ■ When a vortex lattice moves 
coherently at average velocity v through a pinning poten- 
tial, a microscopically periodic velocity modulation |?]] 
is induced at the washboard frequency fi nt = v/a and, 
generally, at integer multiples qfmt- Here, a is the lat- 
tice period in the direction of motion. In presence of an 
rf-force with frequency / these modulations get mode- 
locked when fint/ f = p/q where p, like q, is an integer, 
giving rise to interference plateaus in the 7- V curves. The 
appearance of this phenomenon was shown to depend 
sensitively on the structure of the moving lattice || . 

In this Letter we report the results of mode- locking ex- 
periments on vortex matter driven through mesoscopic 
channels. The channel device j9j consists of an amor- 



phous (a-JNhzGei-z film (x ~ 0.3, T c = 2.68 K and 
thickness d = 550 nm) with a NbN film (T c = 12 K, 
d = 50 nm) on top. Straight parallel channels (width 
W = 230 nm and length 300 fim) were etched in the NbN 
layer. Pinning in a-NbGe is very weak, while in NbN it 
is very strong j[o|, providing easy vortex flow channels 
with strong pinning channel edges (CE's) (see lower inset 
to Fig. Q). Importantly, the structure of vortex matter 
in the channel can be tuned on the 'atomic' scale by 
changing, through the applied field, the (mis)matching 
conditions between the lattice period and channel width. 

In our system the dc-critical current at which vor- 
tices in the channel start to move is determined by 
the (shear) interaction with pinned vortices in the CE's. 
Phenomenologically, J c follows from F c = J C B = 
2Acqq /w e f f |9(] , with F c the critical force density, B the 
induction, cqq the vortex lattice shear modulus and w e f f 
the effective channel width. The parameter A depends 
crucially on the microscopic structure of the vortex ar- 
rays. In fields for which the arrays in the channel and 
the CE's are commensurate (a = ao and w e ff/bo integer 
in the inset to Fig. 0with a = 2& /V3 = (2$ /v / 3B) 1 / 2 
the equilibrium lattice spacing), one expects a large flow 
stress, whereas it should be small for incommensurate 
fields due to the occurrence of misfit dislocations | pd| . 
This picture seems correct, since J c shows oscillations 
and is globally proportional to c 66 However, recent 
studies [ p"2| indicate that this picture is very sensitive 
to positional disorder of the vortex arrays in the CE's, 
particularly when this disorder is sufficient to cause spon- 
taneous formation of dislocations. Recent STM experi- 
ments [ p"3"[ on vortices in NbN indeed show strong devia- 
tions from the regular lattice arrangement. 

This new ingredient turns out to have intruiging con- 
sequences. By means of the mode-locking technique we 
were for the first time able to determine unambiguously 
the number of moving vortex rows, the evolution of the 
microstructure and, from simulations, the microscopic 
origin of the flow-impedance oscillations in the channels. 
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The measurements were performed in a conventional 
four-probe geometry. The rf-current was applied through 
an rf-transformer with balanced transmission lines and a 
matching circuit. We used frequencies from 1 to 200 MHz 
and amplitudes I r f < 5 mA. To avoid heating, the sam- 
ple was immersed in superfluid 4 He. All data shown here 
was obtained at T = 1.9 K. Since strong field- history ef- 
fects exist ||, we focus for simplicity on results of field 
down measurements where first a field ~ /io-ff C 2 of NbGe 
is applied after which we slowly (ngdH/dt ~ —2 mT/s) 
sweep to the field under consideration. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 



I (mA) 

FIG. 1. Main panel: dc-IV curves at 60 mT measured with 
superimposed 6 MHz rf currents (amplitude 4.7, 3.8, 3.0, 2.4, 
1.9, 1.5, 1.2, 0.94, 0.75, 0.59, 0.47, 0.38, 0.30, 0.24, 0.19, 0.15 
and 0mA from left to right). Interference voltages V p are 
indicated. Upper inset: frequency dependence of the lowest 
two interference voltages. Lower inset: schematic geometry 
of the sample and a single channel. Strong pinning CE's are 
marked gray. The (local) vortex lattice parameters and (ef- 
fective) channel width are indicated. 

Figure |l| shows a series of TV curves measured at 60 
mT with superimposed 6 MHz rf-currents of various am- 
plitudes. The dc-/y curve (I r f — 0) shows a nonlinear 
voltage upturn from which I c was determined using a 
100 nV criterion. The application of sufficiently large rf- 
currents leads to clear mode-locking steps at equidistant 
voltage levels V p which are indicated for p — 1,2,3. As 
shown below, q = 1 and thus p — 1 refers to the funda- 
mental and p = 2, 3 to higher harmonics. These voltage 
levels depend linearly on /, as follows for p = 1, 2 from 
the upper inset to Fig. |l| The interference step width AI 
for each harmonic shows Bessel function-like oscillations 
with I r f which will be discussed in detail elsewhere. 

Next we turn to the field dependence of the mode- 
locking phenomenon. The characteristics are better dis- 
played by plotting the differential conductance dl/dV 
versus V, as done in Fig. || for the / = 6 MHz data. 
The results shown are obtained at constant I r f ~ 10/ c in 



fields from 300 mT (bottom) to 60 mT(top) with steps of 
10 mT. The interference steps, now appearing as peaks 
in dl j dV , change in amplitude and in voltage position 
with field. Focusing on the fundamental peak at V\ , it is 
seen that in certain field intervals V\ is field independent, 
e.g. from 60 - 100 mT V x is 4.6 /iV. At 90 mT another 
fundamental peak appears at a voltage of 6.9 /iV, which 
lies between the p — 1 and p = 2 interference peaks in 
lower fields. The second harmonic now coincides with 
the p = 3 peak at 60 mT. Between 100 and 110 mT, the 
lower voltage p = 1 peak vanishes while the 6.9 /xV peak 
(and its higher harmonics) grow and remain detectable 
up to w 180 mT. 
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FIG. 2. Differential conductance dl/dV versus V as ob- 
tained from the IV curves (as in Fig. |y): the applied field is 
varied from 300 mT (bottom) to 60 mT (top) in steps of 10 
mT. The amplitude of the superimposed 6 MHz rf-current is 

Irf/Ic ~ 10. 

The described sequence of appearances and disappear- 
ances of the peaks continues for larger fields. The re- 
sulting variation of Vx with field is shown in more detail 
in Fig. ||(a). We plot the ratio of V\ over / to also 
include data taken at / = 60 MHz. Both data sets col- 
lapse on distinct voltage plateaus separated by abrupt 
steps. Around the steps two consecutive fundamental in- 
terferences can coexist. This behavior is independent of 
frequency and I r f. The voltage separation between ad- 
jacent plateaus corresponds to the value of V\/f of the 
lowest plateau. Therefore we labelled the plateaus by in- 
tegers n. This staircase-like behavior suggests that Vx/ f 
reflects a process in which the number of vortex rows in 
each channel changes with field in a discrete manner. 
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FIG. 3. (a) Vi/f in units $o versus magnetic field /j,qH for 
6 MHz (•) and 60 MHz (■). (b) The normalized width of the 
interference current steps AI/I C of the Vi peaks at 6 MHz. 
Data belonging to odd n are cross-marked for clarity. A7 is 
defined by the gray area in the inset, (c) The field depen- 
dence of the depinning current I c determined from a 100 nV 
criterion. 

The relation between interference voltage and number 
of moving rows can be derived as follows: the general 
interference condition for the velocity of an ordered array 
with average longitudinal vortex spacing a and row spac- 
ing b (see inset to Fig. |l|) is v Piq — (p/q)fa. Thus the 
interference voltage per channel V p ^ q = w e ffBv Ptq , with 
B = $o/a6 the induction, is: 



V p , q = -®ofw e ff/b = ^ofn, 



(1) 



where w e ff/b = n is the number of moving rows. Taking 
into account that we measure the voltage V over ~ 200 
channels, the height of the jumps in Vi/f in Fig. ||(a) cor- 
responds to $o, validating the identification of the peaks 
at Vi as fundamental. We note that, regardless of field, 
no subharmonic peaks (q ^ 1) could be resolved. Such 
subharmonics have been observed in dc-rf driven charge 
density waves and bulk vortex lattices Here, their 

absence might be explained by the fact that our pinning 
potential, which stems from the interaction with pinned 
vortices in the CE's with average spacing ~ et , predom- 
inantly contains modes with wave vector < 2n/a 

The data in Fig. ||(a) combined with Eq. (Q) shows 
unambiguously that in those parts of the flow where vor- 
tices show mode locking, the number of moving rows n 



changes with field by one. We consider now in more de- 
tail the field dependence of the amplitude AI of the mode 
locked steps, i.e. the area under the interference peaks 
(see the inset to Fig. ||(b)). Figure |(b) shows AI as- 
sociated with the voltage plateaus for consecutive n in 
(a) , normalized by I c . At fields around the middle of the 
plateaus, AI/I C is large. This indicates a large degree 
of coherency in the motion and implies that a unique, 
integer number of rows n is moving, corresponding to a 
matching configuration as in the lower inset to Fig. [I]. 
Approaching the edges of the plateaus, AI/I C decreases, 
signalling a reduction of the coherency in the moving ar- 
rays. At the plateau edges, AI/I C is strongly suppressed, 
even if I r f is varied. The coexistence of a residual signal 
from n and nil row configurations at the transition 
fields evidences that these fields correspond to the max- 
imum mismatch situation. 

It is clear that upon approaching a mismatch condi- 
tion, a moving n row configuration experiences increas- 
ing lattice strains which will progressively induce misfit 
dislocations in the structure. Let us now compare the 
observed transitions to the oscillations in I c in Fig. ||(c). 
It is seen that maxima in I c appear at the n — > n ± 1 
transition fields, while the minima in I c are located near 
matching fields ]l5j . Therefore we conclude that, upon 
increasing the density of misfit dislocations in the chan- 
nel, I c is enhanced. This result contrasts the picture in 
which dislocations lower the flow stress and is strongly 
reminiscent of the mechanism expected for the traditional 
peak effect in a random potential. 

As mentioned, the positional disorder of vortex arrays 
in the CE's, which is quenched for fixed field, plays a 
crucial role for understanding the relation between mis- 
match and the dc-depinning current I c . We address this 
issue now in more detail by molecular dynamics simula- 
tions (ll]jlj] of dc-driven channel vortices in the geome- 
try of the lower inset to Fig. ^. Vortex interactions were 
modelled by the London potential with X/ao = 2 (A is 
the penetration depth) The CE vortices were as- 

signed random shifts d with respect to the ideal lattice 
configuration, such that y ((V • d) 2 ) ~ 0.12. We keep 
the average orientation of the CE arrays with a princi- 
pal axis parallel to the channel. The vortex density in 
the channel and the CE's were both (aofro) -1 , while the 
matching condition was tuned by the ratio w e ff /bo. 

First we discuss the results of simulations around 
matching of a 3 row configuration {w e ff /bo = 3.05). The 
(rescaled) J c obtained in such simulations attains a value 
J c = 2Ac 66 /(Bw eff ) = {A/A°)J° with A ~ 0.2,4° fr| 
where A = is the theoretical value for shear 

flow of perfectly matching vortex arrays along ordered 
CE's To understand the origin of this reduced flow 
stress we turn to Fig. ^(a) and (b). In (a) we show 
vortex trajectories during motion over 5ao for a dc-drive 
J ~ 3J°. The trajectories are essentially parallel to the 
CE's and clearly show motion of 3 rows. A Delauney 
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triangulation (Fig. [|(b)) °f a snapshot of (a), shows that 
vortices in the channel are mostly sixfold coordinated 
and elastically coupled. However, at the CE's, disloca- 
tions are present with Burgers vectors (denoted by ar- 
rows) along the flow direction. These dislocations, which 
are mainly induced by disorder in the CE's, are respon- 
sible for the reduction of J c with respect to J®. 



(a). 





ulations is that on approaching a matching state, the 
fault zones gradually heal, yielding a gradual decrease 
of the critical current. Simultaneously the size of coher- 
ent n-row regions grows, providing an explanation for 
the experimentally observed increase in the rf-dc inter- 
ference signal. On changing the field, vortex dynamics in 
the channels thus exhibits a series of smooth structural 
transitions from quasi- ID coherent motion to quasi-2D 
disordered flow in a disordered potential. 

In summary, we investigated the flow of vortices in 
mesoscopic channels bounded by pinned vortices. The 
rf-dc interference measurements yield unambiguous in- 
formation on the field evolution of the number of mov- 
ing rows and coherency in the channel. The depin- 
ning current oscillations exhibit maxima at mismatch 
conditions. Simulations show that this behavior origi- 
nates from blocking of misoriented dislocations in quasi- 
static fault zones, which gradually heal on approaching 
a matching state. 

This work was supported by the 'Stiching voor Funda- 
menteel Onderzoek der Materie' (FOM). N.K. was sup- 
ported by JSPS. V.M.V. was supported by the DOE Of- 
fice of Science. 



FIG. 4. (a) Flow trajectories during motion over 5ao 
around matching (w e ff/bo = 3.05) for J ~ 3Jc. Note the 
positional disorder in the CE arrays, (b) Delauney triangula- 
tion of one snapshot of (a). Direction of the drive / = J x $o 
is shown and Burgers vectors of the dislocations are indicated 
by arrows, (c) Trajectories at mismatch (w e ff/b a — 3.55) 
for J ~ 3.5Jc- (d) Triangulation of a snapshot of (c). All 
data were taken for a channel of length L = 150ao and cyclic 
boundary conditions. 

To study a mismatch situation we consider the case 
w e ff/bQ — 3.55. The critical current we find is character- 
ized by A — 0.065 ~ 0AA°, indeed enhanced compared 
to matching. Figure ^(c) shows the corresponding tra- 
jectories for J ~ 3.5J°. Here, one observes coexistence 
of 3 and 4 moving rows and in addition substantial trans- 
verse displacements and switching of vortices between 
rows. Importantly, the 3 and 4 row regions remain quasi- 
static during motion and their location is determined by 
static disorder in the CE's. The triangulation in (d) of 
a snapshot of the flow reveals numerous dislocations in- 
side the channel. In contrast to the matching case, their 
Burgers vectors assume all possible orientations and are 
mostly misoriented to the flow direction. The motion 
of these misoriented dislocations is blocked by the in- 
teraction with disorder in the CE's and with neighbor- 
ing misoriented dislocations. This results in quasi-static 
fault zones where during the flow the dislocation pat- 
tern changes irregularly and the vortex trajectories are 
jammed. These phenomena form the mechanism under- 
lying the enhanced flow impedance at mismatch. 

The overall picture that has emerged from these sim- 
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